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ABSTRACT 

We show that in order to minimize the uncertainties in the N and O abun- 
dances of low mass, low metallicity (O/H < solar/5) emission-line galaxies, 
it is necessary to employ separate parameterizations for inferring T e (N + ) and 
T e (0 + ) from T e (0 +2 ). In addition, we show that for the above systems, the ion- 
ization correction factor (ICF) for obtaining N/O from N + /0 + , where the latter 
is derived from optical emission-line flux ratios, is (ICF) = 1.08 ± 0.09. These 
findings are based on state-of-the-art single-star H II region simulations, employ- 
ing our own modeled stellar spectra as input. Our models offer the advantage of 
having matching stellar and nebular abundances. In addition, they have O/H as 
low as solar/50 (lower than any past work), as well as log(N/0) and log(C/0) 
fixed at characteristic values of -1.46 and -0.7, respectively. The above results 
were used to re-derive N and O abundances for a sample of 68 systems with 12 + 
log(0/H) < 8.1, whose de-reddened emission-line strengths were collected from 
the literature. The analysis of the log(N/0) versus 12 + log(0/H) diagram of 
the above systems shows the following: (1) the largest group of objects forms the 
well-known N/O plateau with a value for the mean (and its statistical error) of 
-1.43 (+.0084/-. 0085); (2) the objects are distributed within a range in log(N/0) 
of -1.54 to -1.27 in Gaussian fashion around the mean with a standard devia- 
tion of a =+.071/-. 084; and (3) a x~ sc L uare analysis suggests that only a small 
amount of the observed scatter in log(N/0) is intrinsic. 



Subject headings: galaxies: abundances — galaxies: evolution — galaxies: irregular— 
HII regions 
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1. INTRODUCTION 

Low metallicity systems (or metal-poor systems) are low mass emission-line galaxies 
with 12 + log(0/H) roughly in the range 7.2-8.1, where O/H is the oxygen to hydrogen 
number density ratio. They are referred to in the literature as dwarf irregular galaxies 
(dls), their bursting derivatives, blue compact galaxies (BCGs), or H II galaxies, a subclass 
of BCGs dominated in the optical by H II region-like spectra (see Kunth & Ostlin 2000). 
We use oxygen as a measure of metallicity in the above systems because it is the most 
abundant metal in the interstellar medium (ISM), and because its abundance can be easily 
determined from the emission-lines present in their optical spectra. Note that if we adopt 12 
+ log(O/H) = 8.66 ± 0.05 for the solar value (Asplund 2004), then low metallicity systems 
have O/H in the range solar /30-solar/4. The lower limit corresponds to the value of I Zw 18, 
which is the most metal-poor dwarf galaxy ever observed. The upper limit ensures that the 
most metal rich objects are still below the onset of the observed N/O increase with O/H. In 
addition, objects with 12 + log(0/H) below 8.1 have sufficient strength in the auroral line 
[O III] A4363, which is necessary in order to determine the electron temperatures required 
to obtain N/O. 

In the log(N/0) versus 12 + log(0/H) diagram (e.g., Fig. 2 of Pilyugin et al. 2004), low 
metallicity systems form what is known as the primary N/O plateau. For these objects, N/O 
is independent of O/H as is expected for the ratio of two primary elements 1 . The analysis 
of the properties of this plateau plays an important role in understanding the process of 
primary nitrogen production, which is essential for determining the stellar mass range most 
responsible for N-production, constraining stellar N- yields (in particular for massive stars), 
obtaining information about the initial mass function (IMF) of galaxies, and perfecting 
galactic chemical evolution models. Note that although one could think of other candidates 
for studying N/O at low metallicities, such as stars in metal-poor globular clusters, Galactic 
halo stars, or damped Lya systems, which are quasar radiation absorbers with O/H ~ 
solar/10, it is much harder to measure N and O abundances in such systems. This explains 
why the bulk of the N/O data are from low metallicity systems. 

There is disagreement among authors on the vertical thickness of the plateau and on 
the origin of its properties. Pagel (1985) called attention to the large scatter in N/O at fixed 
O/H. However, Thuan et al. (1995), Izotov et al. (1999), and Izotov et al. (2001), found 
that the plateau is narrow, with a dispersion of only ± 0.02 dex in log(N/0), for objects 
with 12 + log(0/H) < 7.6. In addition, the latter authors explained the plateau in terms of 



1 An element is referred to as primary if its production process is independent of the initial metallicity of 
the progenitor star. 
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primary nitrogen from rapidly evolving massive stars (M > 8 M Q ). In contrast, Henry et al. 
(2000) and Chiappini et al. (2003) concluded that these objects are dominated by primary 
nitrogen from intermediate mass stars, which have masses in the range 4-8 M and evolve 
much slower. 

In order to correctly interpret the plateau properties, it is important to up-date the tech- 
niques employed to determine the N/O values. This can be accomplished through modeling. 
In particular, simulations are required to obtain N/O for two reasons. 

Firstly, they are necessary for finding parameterizations for the electron temperatures of 
the regions where N + and + are emitting, i.e., T e (N + ) and T e (0 + ). Unfortunately, because 
the strengths of the auroral features [N II] A5755 and [O II] AA7320, 7230 are commonly 
unavailable in observational data, the above temperatures cannot be obtained from emission- 
line flux ratios of the form, nebular line(s) / auroral line(s), although they are required in 
order to determine the abundances of the above two ions. Since photoionization models 
can provide the auroral lines missing from observations, these lines can be used in turn to 
compute T e (N + ), T e (0 + ), and T e (0 +2 ), independent of each other. Because objects with 
oxygen abundances below 8.1 usually have sufficient strength in the auroral line [O III] 
A4363, T e (0 +2 ) can be directly determined from the observed temperature sensitive line flux 
ratio (^4959 + /5oo7)/-^4363 (Shields et al. 1981). In general, a parameterization of the form 
T e (N + ) = T e (0 + ) = /[T e (0 +2 )], inferred from models, is used to compute N/O at low O/H 
(e.g., Kobulnicky & Skillman 1996; Izotov et al. 1999; Melbourne et al. 2004). Examples of 
this type of relation, based on photoionization models by Stasihska (1990), can be found in 
Pagel et al. (1992) or Izotov et al. (1994). A more recent parameterization based on models 
by Stasihska & Leitherer (1996) is used in Stasihska & Izotov (2003), although the stellar 
spectra used as input to these models reflect the state-of-the-art in 1995. In addition, the 
lowest stellar atmosphere metallicity that they used is solar/10. Ever since, more complete 
synthetic stellar spectra have been published for metallicities as low as solar/20 (e.g., Smith 
et al. 2002). Parameterizations of the form T e (N+) = / 1 [T e (0+ 2 )] and T e (0+) = / 2 [T e (0+ 2 )] 
for systems with metallicities as low as that of I Zw 18 should be obtained and compared to 
each other, using results from up-to-date stellar and nebular models. 

Secondly, models are required for finding the ionization correction factor (ICE) necessary 
for obtaining N/O from N + /0 + . Based on simulations described in Garnett & Shields 
(1987), which make use of Mihalas (1972) and Kurucz (1975, 1979) stellar atmospheres with 
effective temperatures ranging from 38 kK to 55 kK, Garnett (1990) showed that 0.8 < 
(N/0)/(N + /0 + ) < 1.0 in ionized nebulae with one-tenth (Anders & Grevesse 1989) solar 
abundances. More recently, Izotov et al. (2004) confirmed that N/O N + /0 + is a reasonable 
approximation in metal-poor H II galaxies. They established this using photoionization 
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models of Stasinska & Izotov (2003) with input stellar spectra generated with the stellar 
atmosphere code CoStar (Schaerer & de Koter 1997), as well as simulations computed with 
the evolutionary synthesis code Starburst99 (Leitherer et al. 1999), using Smith et al. (2002) 
stellar spectra generated with WM-Basic (Pauldrach et al. 2001). Note that the spectra 
of Smith et al. (2002) are available for metallicities down to Z Q /20. Unfortunately, the 
uncertainty in the ICF is not given in Izotov et al. (2004), however, the value of the ICF is 
sensitive to the ionization parameter 2 . This parameter is a poorly constrained variable (e.g., 
log(U) ranges from -5 to -1 in models by Garnett 1989). Note that more recently, Stasinska 
& Leitherer (1996) found that only a small range for U can fit the emission-line diagnostics 
from the few H II galaxies that they modeled. However, as these authors conclude, U is 
yet to be observationally constrained. In addition, the value of the ICF also depends on the 
method employed to compute N + /0 + . In particular, if the latter is derived from optical 
line flux ratios, the ICF will depend on the temperatures used to derive the N + /H + and 
0+/H+ ratios. 

Our goal is to re-derive N and O abundances for a carefully selected sample of 68 low- 
metallicity systems whose de-reddened optical emission-line strengths were taken from the 
literature (Campbell et al. 1986; Walsh & Roy 1989; Pagel et al. 1992; Skillman & Kennicutt 
1993; Izotov et al. 1994; Thuan et al. 1995; Kobulnicky & Skillman 1996; Izotov et al. 
1997; Kobulnicky & Skillman 1998; Izotov & Thuan 1998a,b; Izotov et al. 1999; Melbourne 
et al. 2004). This was accomplished using independent temperature parameterizations for 
T e (N + ) and T e (0 + ) (which diverge at low metallicities, see below) and applying an ICF 
for obtaining N/O from N + /0 + . Our temperature patameterizations and ICF are based 
on results from a representative grid of photoionization models computed with the code 
CLOUDY version 96.01 (Ferland et al. 1998). These employ state-of-the-art modeled stellar 
spectra that we generated with PHOENIX version 13.08.04A (Hauschildt & Baron 1999, 
2004; Aufdenberg et al. 2002). The purpose of our work is to carefully study the morphology 
of the N/O plateau, assess the uncertainites, and then quantify the amount of natural scatter 
among plateau objects. 

We discuss our general procedure for deriving the number density ratios O/H and N/O 
from optical emission-lines in § 2.1. In § 2.2 we describe our models and explain how we ob- 
tained the electron temperatures, electron density, and ICF relevant to the above abundance 
determinations. Our method for determining abundance uncertainties is explained in § 3, in 
which we also compare our abundances with literature results. In section § 4 we analyze the 



2 The ionization parameter can be denned as U = Q(H) / [ 4 7r i 2 n(H) c ], where, r Q is the distance from 
the ionizing source to the illuminated face of the cloud, n(H) is the total hydrogen density, c is the speed of 
light, and Q(H) is the stellar emission rate of hydrogen ionizing photons (Ferland et al. 1998). 
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scatter in the N/0 plateau. Finally, section 5 gives a summary of our work. In a companion 
paper to this one (Henry et al. 2006), chemical evolution models are combined with Monte 
Carlo techniques to further interpret the plateau morphology. 

2. METHOD 

2.1. General Procedure 

The N and O ionic abundances for our sample of low metallicity systems were determined 
directly from published optical emission-lines using the general expression: 

H+ I Hf) X u (Te,N e ) At 4861' 1 j 

where the fraction on the left is the number density ratio of the N or O ion relative to H + , the 
first fraction on the right is the line flux ratio of a nebular forbidden feature of ion X* relative 
to the strength of N e is the electron density (cm -3 ), T e is the electron temperature 
of the region where the relevant ion is emitting (K), a e J^ is the effective recombination 
coefficient of H/3 which includes radiative and three-body processes (cm 3 s _1 ), \ u is the 
fraction of ions X* with an electron in the upper level of the transition of interest, A" 
is the corresponding spontaneous de-excitation rate coefficient (s _1 ), and the last term is 
the wavelength ratio of the line of interest and H/3 (A). Note that equation (1) is based 
on the assumption that H(3 arises from recombination. The contribution to H lines from 
collisional excitation was neglected because the excitation potentials of H levels are much 
higher than the average thermal equilibrium temperature that characterizes H II regions 
(Osterbrock, D. E. & Ferland, G. J. 2006). In addition, equation (1) assumes that nebular 
forbidden lines originate from collisionally excited levels. Collisional excitation is significant 
in this case because the low-lying energy levels of the relevant ions are of the order of kT. 
However, according to Rubin (1986), recombinations of +2 can excite the nebular doublet 
[O II] AA3726,3729, used to compute + /H + . The effect of this process will be analyzed in 
the future. 

Although the presence of the nebular He II 4686 A emission-line in the spectra of 
several metal-poor systems implies the presence of unobserved +3 /H + , the contribution 
of this ion to the total oxygen abundance amounts to only a few percent according to our 
photoionization models. Therefore, we obtained the oxygen abundance by assuming that 
O/H = ( 0+ + 0+ 2 ) / H+. The N/O ratio was determined from N/O = ICF x N+/0+, 
where the ICF was derived using models as explained in § 2.2. 
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A 5 th order fit to data from Storey, P. J. and Hummer, D. G. (1995) was adopted for 
a^p, i.e., = where £ = Yli=o a (^e ■ The parameters are: a = —12.404592, 
a x = -3.47193796 x 10" 4 , a 2 = 4.98365006 x 10~ 8 , a 3 = -3.77545451 x 10~ 12 , a 4 = 
1.33944026 x 10~ 16 , and a 5 = -1.75120267 x 10~ 21 . The fraction \ u was determined using 
a 5-level atom routine written and tested thoroughly for accuracy by one of us (A. N.). The 
equilibrium equations for each electronic level were solved using Gauss- Jordan elimination 
and the maximum pivot technique (Chapman 1998). Column (1) of Table 1 shows the list 
of ions involved in the determination of N/O. Columns (2), (3), and (4) of the same table 
give references for the energy levels, transition rates, and collisional strengths, respectively, 
employed in the 5-level atom calculations. Note that S + is included in Table 1 since N e ([S 
II]) was employed in our calculations as explained in the next section. 



2.2. Electron Temperatures, Density, and ICF 

2.2.1. Electron Temperatures 

In equation (1), three parameters, a e M, x"? an d N e , depend on T e . In theory, T e (N + ), 
T e (0 + ), and T e (0 +2 ) can be obtained from the T e sensitive line flux ratios, R[N II] = 

(-^6548 + -^6583) /-^5755, -R[0 II]= (-Z3726 + -^3729)/(-^7320 + ^723o) , and _R[0 III] = (^4959 + -^5007) /-^4363, 

respectively (Shields et al. 1981). Unfortunately, in practice T e (N + ) and T e (0 + ) cannot be 
directly inferred from their corresponding ratios because the auroral line [N II] A5755 and 
the auroral doublet [O II] AA7320, 7230 are usually absent in observational data. However, 
photoionization models can predict the auroral line strengths required for deriving theoreti- 
cal relations of the form T e (N+) = ,A[T e (0 +2 )] and T e (0+) = / 2 [T e (0+ 2 )]. The independent 
variable of these relations is T e (0 +2 ), since this temperature can be inferred directly from the 
observed R[0 III] ratio of low metallicity systems as long as [O III] A4363 is detected with 
good S/N. Note that the advantage of finding individual parameterizations for T e (N + ) and 
T e (0 + ) is that this minimizes the uncertainty in N + /0 + (and N/O) which is introduced by 
assuming that T e (N + ) and T e (0 + ) are equivalent functions of T e (0 +2 ), as discussed in § 1. 

In order to find separate temperature parameterizations for T e (N + ) and T e (0 + ), deter- 
mine the ICF for obtaining N/O from N + /0 + , and carefully estimate uncertainties in N/O, 
we computed a grid of photoionization models, representative of low metallicity systems. In 
the following paragraphs, we introduce PHOENIX, the stellar atmosphere code used to gen- 
erate the spectra that served as input to our CLOUDY photoionization models, we explain 
the choice of parameter space covered by our stellar and nebular simulations, we describe 
how we proceeded to obtain the temperature parameterizations and ICF from our models, 
and we derive the electron temperatures, electron densities, and ionic abundances of our 
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sample objects. 

PHOENIX is a general stellar atmosphere code which operates in one-dimensional spher- 
ically symmetric or plane-parallel geometry. It can handle static atmospheres (e.g. stars, 
Jupiter-like planets, brown dwarfs) as well as moving atmospheres (supernovae, novae, mas- 
sive hot stars, etc.). The code simultaneously solves the equations of radiative transfer (in 
the proper geometry), structure, and level population in Non-Local Thermodynamic Equi- 
librium (NLTE). 

We used PHOENIX's wind module package (Aufdenberg et al. 2002) to generate syn- 
thetic spectra in the range from 100 to 8xl0 6 A, for six windy O-dwarf stars with solar/5, 
solar/20, or solar/50 metal abundances, except for N and C, which we scaled in order to be 
consistent with typical proportions found in low metallicity systems, i.e., log(N/0) = -1.46 
and log(C/0) = -0.7 (Henry et al. 2000). The set of solar abundances was modified in order 
to match that of the code we used for our photoionization models 3 . Other input parameters 
to the models were based on those employed by Smith et al. (2002), who generated WM- 
Basic synthetic spectra for the same type of stars, with solar/5 and solar/20 proportions for 
all metals (including N and C). 

We adopted a /3-law for the wind velocity, i.e., 

v(r) = Voo ( 1 - R*/r f, (2) 

where r is the distance from the center of the star, is the radius of the "photosphere" 
(considered to be the base of the wind), and is the terminal velocity, assumed to be 
reached at 100 x R*. The /5-index was assumed to be 0.8 (Lamers & Cassinelli 1999). 
The density structure in the dynamic region (v(r) > 0) was calculated from the continuity 
equation: 

M , , 

p(r) = - — --, (3) 

4 Ti r 2 v(r) 

where, p{r) is the mass density and M is the mass loss rate. In the hydrostatic region 
(v(r) = 0) the density structure was calculated assuming hydrostatic equilibrium. All models 
have 45 zones for the dynamic region and five zones for the hydrostatic region. 

A great deal of acceleration comes from resonance line radiation pressure of metals. 
Therefore, as in Smith et al. (2002), the dependences of and M on the metallicity were 
parametrized by: 

^oo = ^oo,0 Z°- 13 , (4) 



3 Our software uses solar abundances by Allende Prieto et al. (2001). 
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and 

M = M Q Z - 8 . (5) 

Values for Voo,e and M Q can be found in Smith et al. (2002). A summary of the parameters 
characterizing our PHOENIX models is given in Table 2. For each model listed in column 
(1), columns (2) through (8) show the following information: (2) atmosphere's metallicity, 
(3) stellar mass loss rate, (4) wind's terminal velocity, (5) photospheric temperature, (6) 
surface gravity, (7) stellar radius, and (8) total luminosity in the range specified within 
parenthesis. In Figure 1, we display the ionizing portion of PHOENIX (solid lines) and 
WM-Basic (dashed lines) spectra for two solar/20 O-dwarfs. The effective temperatures of 
the stars are 43 kK (top panel) and 50 kK (bottom panel). Note that the PHOENIX models 
differ from the WM-Basic models in that the former have more hydrogen ionizing photons 
and fewer helium ionizing photons. This is due to the fact that PHOENIX absorbs the 
He ionizing flux due to the presence of more lines. Therefore, the normalized luminosity is 
re-distributed towards H ionizing energies. 

We used CLOUDY to simulate the spectra of 24 single-star H II regions. The PHOENIX 
star used as input to each model is listed in column (1) of Table 3. Note that we reduced the 
spectral resolution of the PHOENIX energy distributions in order to fit CLOUDY's default 
resolution. This was done using the CONGRID routine implemented in IDL version 5.5 4 . 

The following general characteristics are shared by all of our models: spherical geometry 
(ionizing source fully covered), constant H density n(H), homogeneous distribution of the 
gas (filling factor of 1), and static configuration. Studying the effects of using more complex 
models falls beyond the scope of this paper. 

All nebulae have He/H = 0.1 and metal abundances identical to those of the atmospheres 
of the illuminating stars, except for the refractory elements Na, Mg, Al, Si, Ca and Fe, which 
are depleted by a factor of 10 since they are likely to be found in the form of grains in H 
II regions (Savage & Sembach 1996). Accordingly, grains were included in our simulations 
by turning on the "grains ISM" feature implemented in CLOUDY, with their abundances 
scaled down by a factor of 1/5, 1/20, or 1/50, in order to be consistent with the oxygen 
depletion of the model 5 . 

In order to study the effect of varying the nebula's electron density, we ran models 



4 This Solaris sunos spare data analysis software was developed by Research Systems Inc. (2001). 

5 Note that CLOUDY's ISM grains are based on an older set of solar abundances than the set used by 
the code to scale the H II region abundances. The latter set is based on the abundances given by Grevesse 
& Sauval (1998), Holweger (2001), and Allcnde Prieto et al. (2001, 2002) 
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for two different values of n{H), 10 or 100 cm" 3 , which are typical values measured in low 
metallicity systems. We also varied the distance from the star to the illuminated face of the 
cloud (r ), adopting either 2 pc, (which is bigger than the radius of a planetary nebula) or 
5 pc (which allowed us to cover a range for the ionization parameter comparable to what 
is found in the literature). The above choices resulted in -3.17 < log(U) < -0.42. The 
parameters characterizing our CLOUDY models are summarized in Table 3, which gives 
O/H in column (2), the H density in column (3), r in column (4), log(U) in column (5) and 
the nebular geometry in column (6). Note that the latter is plane-parallel if Ar/r D < 0.1, a 
thick shell if 0.1 < Ar/r Q < 3, and spherical if Ar/r Q > 3, where Ar is the thickness of the 
cloud set by the stopping temperature, chosen to be 100 K. 

In order to find parameterizations of T e (N + ) and T e (0 + ) with respect to T e (0 +2 ), from 
models, T e (N + ), T e (0 + ), and T e (0 +2 ) need to be determined independent of each other first. 
This can be accomplished using two different methods. The first method, corresponding to 
common practice (Garnett 1990; Stasihska 1990; Stasihska & Leitherer 1996; Stasihska & 
Izotov 2003), is to compute each temperature from: 

jT e {A + i)N e N{A + J)dV 



(T e (A +J )) = 



N e N(A+i)dV ' ^ 

i.e., each temperature is an ion and N e weighted volume T e mean. However, the above three 
temperatures can also be obtained from their respective line flux ratios, R[N II], R[0 II], 
or R[0 III], as explained below. This second method has the advantage of being consistent 
with the observational way of obtaining T e (0 +2 ), i.e., from line flux ratios. Ideally, tem- 
perature parameterizations should be independent of the method employed to compute the 
individual temperatures. We will compare temperatures computed with different methods 
further on, but first, we will explain how we determined T e (N + ), T e (0 + ), and T e (0 +2 ) from 
their respective line flux ratios. 

In order to accomplish the latter, we used an iterative procedure based on Brent's 
method for finding roots (Press et al. 1996). For each of the 24 photoionization models, the 
idea is to use our 5-level atom routine to compute repeatedly the theoretical line flux ratio, 
Rthe, for assumed values of T e over a range of 4-25 kK until it matches the model-predicted 
value Rmod to within a tolerance of 0.001. We carried out this procedure separately for the 
two density regimes of 10 and 100 cm" 3 . Using a maximum of 100 iterations to solve the 
equation numerically, we searched for the value of T e which satisfied the relation: 

Rthe[T e ,N e ] — R mo d 



Rthe[T e , N e 



= 0.001. (7) 



Henceforth, we will refer to models with n(H) = 10 cm 3 and n(H) = 100 cm 3 as the 
nlO and nlOO models, respectively. Our temperature parameterizations for the nlO models 
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are 

T e (N + ) = -5950 + 2.256 x T e (0+ 2 ) - 6.28 x 10" 5 x T e (0 +2 ) 2 (8) 

and 

T e (0 + ) = -5470 + 2.131 x T e (0 +2 ) - 5.54 x 10~ 5 x T e (0 +2 ) 2 , (9) 
while for the nlOO models they are 

T e (N + ) = -7720 + 2.438 x T e (0 +2 ) - 6.56 x 10~ 5 x T e (0 +2 ) 2 , (10) 

and 

T e {0 + ) = -7330 + 2.325 x T e (0 +2 ) - 5.82 x 10~ 5 x T e (0 +2 ) 2 . (11) 

For each model, Figure 2 shows two independent values of T e (0 +2 ) on left panel, and 
two independent values of T e (0 + ) on the right panel. On the left panel, the value obtained 
using R[0 III] is given vertically, while the value (T e (0 +2 )) predicted by equation (6) is 
given horizontally. The agreement between the two T e (0 +2 ) temperatures is excellent. On 
the right panel, the parameterized T e (0 + ) given by either equation (9) or (11) [depending 
on the input value of n(H)] is given vertically, with open circles representing nlO models 
and filled circles representing nlOO models. On the same panel, the horizontal value is 
that predicted by equation (6). The agreement between the two T e (0 + ) temperatures is 
reasonable, with the parameterized value being higher than (T e (0 + )) for most models. The 
biggest difference found is T e (0 + ) - (T e (0 + )) of about 800 K. The effect of this difference 
on the N + /0 + ratios is discussed below. Since unfortunately CLOUDY does not print out 
{T e (N + )), the latter temperature was not compared to the parameterized T e (N + ). 

Figure 3 features T e (N+) vs. T e (0+ 2 ) (top panel) and T e (0+) vs. T e (0+ 2 ) (bottom 
panel), where results for the nlO and nlOO models are shown with open and filled circles, 
respectively. In this figure, the solid lines are best fits to models sharing the same value of 
n(H) while the dotted lines are extrapolations to low metallicities of parameterizations of 
the form T e (N+) = T e (0+) = /[T e (0+ 2 )], from Pagel et al. (1992), Izotov et al. (1994), and 
Stasihska & Leitherer (1996). Note that the latter three fits are based on temperatures given 
by equation (6). Henceforth we will refer to them as pag92, izo94, and sta96. 

The following can be noticed when analyzing Figure 3. For a given ion (N + or + ), 
T e (vertical axis) at a given T e (0 +2 ) is in general higher for the nlOO models than for the 
nlO models. This is due to the fact that as n(H) increases, U decreases, bringing the regions 
where N + and + are emitting closer to the ionizing source, making them hotter. This 
behavior is illustrated in Figure 4, which has three panels. The top panel shows T e as a 
function of depth from the illuminated face of the cloud for two models with T e ff = 43 kK, 
r Q = 2 pc, and O/H = solar/50, only differing in their value of n(H). The lower two panels 
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show the ionization fractions N + /N (solid line) and + /0 (dashed line) as functions of depth 
for the two models. The middle panel corresponds to the nlO model, while the bottom panel 
corresponds to the nlOO model. The dotted lines clearly show that T e at the peak of N + /N is 
higher by roughly 1000 K for the nlOO model. The same is true for T e at the peak of + /0. 

Also in Figure 3 we see that for a given value of n(H), T e (N + ) is in general lower than 
T e (0 + ) [except for the two models with T e ff = 43 kK, r Q = 2 pc, and O/H = solar/5, for 
which they are almost identical]. This is shown more clearly in Figure 5, which features 
T e (N + ) vs. T e (0 + ) for all models, and is a result of the fact that the regions where N + and 
+ emit do not completely overlap, the + region extending more towards the star into 
hotter gas, as we pointed out earlier in our discussion of Figure 4. Our models also indicate 
that this offset becomes larger with rising stellar temperature and falling metallicity. 

In order to show that it is better to use T e (N + ) independent of T e (0 + ) when determining 
the N + /0 + ratios of low metallicity systems, we derived the N + /0 + ratios of our sample 
objects using three different methods but always adopting N e = 100 cm -3 . Methods 1 and 
2 both assume that T e (N + ) = T e (0 + ), the assumption frequently made by researchers in the 
past. In method 1, both temperatures are computed using equation (10), while in method 2 
equation (11) is used instead. Then values of N + /H + , + /H + , and ultimately N + /0 + are 
determined from these derived temperatures. Finally, method 3 uses equations (10) and 
(11) to compute N + /H + , + /H + , respectively. The resulting three sets of N + /0 + ratios are 
shown in Figure 6. The top panel is a plot of method 1 versus method 3, while the bottom 
panel is a plot of method 2 versus method 3. We see that using either method 1 or 2 results 
in systematically underestimating the N + /0 + ratios by roughly .05-. 10 dex. Therefore, in 
order to minimize uncertainties in determining N + /0 + in low metallicity systems, one clearly 
should derive T e (N+) and T e (0+) from T e (0+ 2 ) separately. 

A comparison of our temperature parameterizations (eqs. [8] and [9]) with those of 
pag92 and izo94 in Figure 3 shows good agreement among the three functions at a density 
of 10 cm -3 , where the pag92 and izo94 formulations are based on photoionization models by 
Stasihska (1990), most of which have n(H) = 10 cm" 3 . On the other hand, for a density of 
100 cm -3 the agreement between our results (eqs. [10] and [11]) and pag92 and izo94 is less 
satisfactory. For the temperature formulation by sta96, there is relatively poor agreement 
with our low density results for T e (N + ) and T e (0 + ), although the agreement is better for 
our results for a density of 100 cm -3 , especially in the case of T e (0 + ). Differences between 
the N + /0 + values obtained with sta96 and our method are due to the fact that sta96 
constitutes an extrapolation to lower metallicites of the published relation. In addition, the 
latter is based on simulations of H II regions created by an evolving group of stars, while 
we use single non-evolving star simulations. Finally, our method uses independent relations 
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for T e (N + ) and T e (0 + )with respect to T e (0 +2 ). To test the impact of the density effect on 
our final N + /0 + values, we used equation (1) to calculate this ratio for each of our sample 
objects for the two density regimes, nlO and nlOO. The results are shown in the upper 
left panel of Figure 7, where the nlOO regime (ordinate) is plotted against the nlO regime 
(abscissa). Clearly, the points closely follow the diagonal line corresponding to equal values 
of N + /0 + in both cases. 

Further comparisons of our temperature method with the three literature results just 
discussed are shown in the remaining panels of Figure 7. In each case we have calculated 
N + /0 + values at a density of 10 cm— 3 for our sample objects using both our temperature 
method (vertical axis) and that of pag92 (upper right), izo94 (lower left), and sta96 (lower 
right), shown on the horizontal axis in each graph. The diagonal lines show the one-to-one 
relation. Note the particularly poor agreement between our method and those of pag92 and 
izo94, especially for N + /0 + values below -1.2. This disagreement can be easily traced to 
the relatively low T e (0 + ) temperatures predicted by pag92 and izo94 when T e (0 +2 ) exceeds 
14000 K, resulting in higher + /H + and lower N + /0 + ratios than ours. On the other hand, 
the apparently closer agreement with sta96 arises because their parameterizations for both 
T e (N + ) and T e (0 + ) are offset from our curves by roughly the same amount in each case, 
resulting in relatively little difference in our computed N + /0 + ratios. 

Since in general the only N e sensitive line flux ratio available from observations is the 
[S II] ratio Iq7\q/Iq7zii we assumed that N e was equal to iV e ([S II]) throughout the observed 
H II regions. For 42 low metallicity systems in our sample the [S II] ratio was detected 
with reasonable S/N. For this group of objects, iV e ([S II]) was computed by applying Brent's 
method as described above, but replacing R mo d with the observed value, Robs- We pro- 
ceeded as follows. We first fixed N e at 100 cm' 3 (a typical value for H II regions) and 
used R obs [0 III] and equation (7) to determine T e (0+ 2 ). Next, from T e (0+ 2 ), R obs [S II] and 
equation (7) we computed N e . This new value of N e was then employed in the derivation of 
the N and O abundances of the above 42 objects. For the remaining objects, we used the 
average N e found for this subsample, i.e., ~200 cm -3 . 

Our final values for N e , T e (0 +2 ), T e (0 + ), and T e (N + ) are listed for each object in our 
sample in columns (2) through (5) in Table 5. 

Finally, we used our models to derive an ICF for obtaining N/O from N + /0 + . This 
was accomplished by computing the average of (N/0)/(N + /0 + ) for the nlOO models, where 
relevant line fluxes and parameterizations were employed to obtain each N + /0 + ratio, while 
N/O was fixed and given by the input value to our models, i.e., log(N/0) = -1.46. The 
average ICF for these models is 1.08 with a standard deviation of 0.09. This contrasts with 
common practice which is to use ICF = 1. 
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Unfortunately, the ICF depends on the method employed to compute N + /0 + . A com- 
mon practice is to compute N + /0 + from the ion fractions (N + /N) and (0 + /0) given by: 



and the input values of N and O. A comparison between the N + /0 + ratios that we used to 
compute our final ICF, and those given by equation (12), is shown in Figure 8. In this figure, 
the difference between log(N + /0 + ) and log((A + /0 + )) is plotted against log((A + /0 + )). The 
latter value is that predicted by equation (12), while log(N + /0 + ) is the value obtained from 
relevant line flux ratios and the temperature parameterizations for the nlOO models. Note 
that Figure 8 illustrates the need of an ICF for obtaining N/O at low metallicities, no matter 
which method is used to compute it. Clearly, log((A + /0 + )) deviates from log(N/0) which 
as mentioned above is fixed at -1.46 for all models. In fact, the use of (N + /0 + ) when 
computing the ICF results in an average of 1.28 with a standard deviation of 0.1. 

The vertical difference in Figure 8 is at most ~0.11 dex. It is positive because N + / + val- 
ues from line flux ratios and parameterized temperatures are higher than those predicted by 
equation (12). This fact is reflected in the smaller ICF obtained when using line flux ratio 
N + /0 + values. Since our method is consistent with the way in which N + /0 + is measured 
observationally, i.e., from line flux ratios, we applied an ICF of 1.08 when deriving the N/O 
abundances of our sample objects. The ICF's contribution to the final N/O uncertainties 
was accordingly taken to be ±0.09. 



We used the above results for density, temperatures, and ICF determinations to derive 
new values for the O/H and N/O ratios for our sample objects. The corresponding uncer- 
tainties were obtained using standard error propagation, where the quantities considered in 
our estimation of the abundance errors is given in Table 4. The quantities listed in column 

(1) were assumed to be functions of the variables given in column (2). Note that because 
N e depends on T e and visa versa, and since in addition, N e had to be assumed for a significant 
number of objects (see below), N e does not appear in Table 4. 

For each object in our sample, Table 5 gives the object name in column (1), N e in column 

(2) , T e (0+ 2 ), T e (0+), and T e (N+) in columns (3), (4), and (5), respectively, and 0+/H+, 
+2 /H + , and N + /H + in columns (6), (7), and (8), respectively. Our final abundance results 
are given in Table 6, which lists object name in column (1), references for the observed 
emission-line strengths in column (2), our derived values for 12 + log(0/H) ± a and log(N/0) 
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± a in columns (3) and (4), and literature values and references in columns (5), (6), and 
(7), respectively. 

Figures 9 and 10 feature comparisons of our O/H and N/O values (vertical axes) against 
published values. The diagonals show points of one-to-one correspondence. Our O/H values 
are systematically lower with respect to past calculations due primarily to our new temper- 
ature scheme for calculating + /H + . However, in general, the error bars show agreement 
between our results and literature values. Our N/O values are also offset (higher by ~ 0.05 
to 0.1 dex) due the use of new temperature parametrizations and an ICF for obtaining N/O. 

In Figure 11, we plot our uncertainties in log(N/0) against corresponding literature 
values. Some of the published uncertainties were estimated using Monte Carlo simulations 
to propagate the errors in the relevant line strengths (Campbell et al. 1986; Kobulnicky 
& Skillman 1996), and thus points occupy both sides of the diagonal due to the random 
nature of the latter technique. We point out that log(N/0) values with small uncertainties 
generally correspond to objects with small uncertainties in relevant emission-line strengths. 
The tendency for our error estimates to systematically exceed published ones is primarily 
due to the contribution of ICF uncertainty for each of our objects. However, we strongly 
argue that inclusion of this additional factor produces more realistic uncertainties than those 
previously published. 



4. ANALYSIS 

Figure 12 is a plot of N/O versus O/H for the abundances and uncertainties listed in 
Table 6. Recall that N/O is a ratio which measures differential nucleosynthesis, while O/H 
is a metallicity gauge and measures the extent of chemical evolution in the area immediately 
surrounding the observed emission-line region. Furthermore, each data point in Figure 12 
provides a snapshot of current conditions within a galaxy, each of which arrived at its present 
position by following a path determined by its own unique star formation history. 

Figure 13 shows the distribution of N/O values for our sample, using a bin size of 
0.02 dex. Two distinct groups of objects are apparent. The larger population comprises 
objects extending over the entire O/H range considered here and possessing N/O values 
which appear to be distributed in near-Gaussian fashion between limits of -1.54 < N/O 
<-1.27. These limits are indicated with dashed lines in Figures 12 and 13. The remaining 
objects with N/O > -1.27 appear to belong to a smaller population which is concentrated at 
the high end of our O/H range and possess relatively high N/O values compared with the 
larger group. 
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The 52 objects making up the larger population form what we will define for purposes 
of analysis as the N/O plateau. The N/0 plateau has been discussed by numerous authors, 
including Garnett (1990); Thuan et al. (1995), Kobulnicky & Skillman (1996), and Pilyugin 
et al. (2003). We argue that by treating this group of objects separately, we are not doing so 
in an arbitrarily fashion. Rather we are led to this result for the following reasons. First, the 
distribution of the entire sample in the N/O dimension shown in Figure 13 strongly hints at 
the existence of two distinct groups of objects. Second, our current theoretical understanding 
of N production strongly suggests that we should expect a narrow, horizontal track at low 
metallcities where the production is primary. Both analytical and numerical models such 
as those presented in Henry et al. (2000) predict the presence of this track. Finally, a 
theoretical analysis of the N/O plateau shape by Henry, Nava, & Prochaska (2006) using 
numerical chemical evolution models finds that an additional N production process triggered 
at 0/H>7.8 is needed to explain the high N/O objects. Thus, we conclude that current 
observational and theoretical evidence suggests that what we have defined as the N/O plateau 
consists of objects which have been subjected to similar N production mechanisms, and that 
we are justified in excluding, for the time being, those objects in the smaller group which 
possess higher N/O ratios. 

We now discuss some of the statistical aspects of our measurements. In doing so we 
note that all means, errors, and standard deviations relating to our sample and given below 
were initially computed in linear space and then converted to their corresponding logarithmic 
quantities. This is not necessarily true for values quoted from other papers, however. 

The log of the weighted mean of N/O for the 52 objects lying on the plateau is -1.43 
(+.0084/-. 0085), where the values in parenthesis represent the error in the mean, not the 
standard deviation, and were derived from observational uncertainties quoted in column (4) 
of Table 6 6 . The value for plateau objects with O/H below 7.8 is -1.44 (+.011/-.012), while 
above 7.8 it is measured to be -1.41 (+.012/-. 013), so there is no evidence of metallicity 
dependence of the mean. These results are in good agreement with those of Izotov et al. 
(1999), who obtained -1.47 (a = +.14) for their total sample, -1.60 (a = +.02) for their 
low metallicity (O/H < 7.6) objects, and -1.46 (+.14) for their high metallicity (O/H > 
7.6) objects. Likewise, our results are in close agreement with Garnett (1990), who found 
log(N/0) = -1.46 (+.10/-. 13) for plateau objects. Finally, the standard deviation about the 
weighted mean for our plateau objects, expressed logarithmically, is +.071/-. 084, nearly an 
order of magnitude larger than the error in the mean. 



6 Note that the error in the mean depends upon the uncertainties in the individual measurements as well 
as the square root of the number of objects in the sample, being inversely related to the latter (see Bevington 
& Robinson 2003). 
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An interesting question pertaining to the distribution of points in the N/O-O/H plane 
concerns the nature of the observed N/O scatter (as measured by the standard deviation) of 
the plateau objects and to what degree it is related to errors in abundance measurements. 
That is to say, is the spread in N/O associated with the plateau due mostly to intrinsic or 
statistical scatter? To investigate this point further, we assumed a Gaussian distribution in 
the parent N/O plateau population, as suggested by the distribution shape in Figure 13, and 
performed a x-square analysis of the observed point distribution of log(N/0) for the plateau 
using our calculated weighted mean value for log(N/0) of -1.43 along with the established 
errors in the abundances, as given in column (4) of Table 6. We found that the value for 
the reduced x-square 7 for the 52 plateau objects is 1.244, consistent with a probability of 
about 90% that our sample of objects could not have been drawn randomly from a parent 
population characterized by our calculated mean log(N/0), error, and standard deviation. 
In other words, some additional scatter is needed to explain the observed spread in N/O 8 . 

In order to estimate the magnitude of this additional scatter, we recalculated reduced 
%-square values for the 52 plateau objects numerous times, each time adding intrinsic scatter 
to the measurement errors in quadrature to simulate real scatter among the objects, with 
a target value for the reduced x _sc l uare of unity. Such a value was reached when a for the 
intrinsic component was 0.024, or about 1/3 of the magnitude of the standard deviation and 
the typical statistical uncertainties in column (4) of Table 6. Therefore, the implication of 
this x- S( l uare analysis is that only a small portion of the scatter observed for the plateau 
objects is intrinsic. Considering that the many of the uncertainties associated with the line 
strengths taken from the literature for our analysis are only estimated and not determined 
in a rigorous fashion, we cannot conclude with any confidence that a significant portion of 
the vertical scatter in N/O is real. Rather, the issue can only be decided by using better 
measurements including well-determined uncertainties in the abundances. 



7 Thc reduced x-squarc, xt — jy-i Si — m] 2 |: where yi and a% are the N/O value and uncertainty 

of the ith object, fi is the sample average, and N is the number of sample objects. Within the summation, 
the denominator represents the measurement errors. 

8 Concerned that a few points with small observational uncertainties could be artificially inflating the 
total x-square value, we calculated each object's contribution to the total x-square and found no points to 
be excessively high. 
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5. CONCLUSIONS 

This paper has been devoted to a detailed review of the method for deriving N and O 
abundances in low mass emission-line galaxies with 12 + log(0/H)< 8.1. These low metallic- 
ity systems are classified in the literature as dwarf irregular galaxies, blue compact galaxies, 
or H II galaxies. The principal motivation for our work was to estimate the contribution of 
the N/O uncertainties to the scatter among N/O plateau objects in the log(N/0) versus 12 
+ log(0/H) diagram. To carry out this study, we selected a sample of 68 objects from the 
literature and used their published de-reddened emission-line strengths to calculate a homo- 
geneous set of abundances, paying special attention to the determination of [N II] and [O II] 
electron temperatures for which the necessary but weak auroral lines are usually absent from 
the data. In such cases these temperatures, which are needed to compute the abundances 
of their associated ions (N + and + ), must be inferred from parametric relations in terms 
of the more accessible [O III] temperature. We established theoretical relations for the rel- 
evant temperatures based on single-star H II region simulations (modeled with CLOUDY) 
and employing our own input stellar spectra (modeled with PHOENIX). These relations 
were derived from line flux ratio temperatures in contrast with common practice, which is 
to use ion and iV e -weighted volume temperature means. We also used our models to infer 
an ICF for obtaining N/O from N + /0 + , and to carefully estimate N/O uncertainites by ap- 
plying standard error propagation methods. This ICF was derived from line flux ratio ionic 
abundances, also in contrast with common practice, which is to use iV e -weighted volume 
ion fraction means. Our method has the advantage of being consistent with the procedure 
employed to derive T e (0 +2 ) and N + /0 + from observations, i.e., from line flux ratios. Note 
that N + /0 + values from line flux ratios and parameterized temperatures are higher than 
those predicted by equation (12). However, this fact is taken care of by our value for the 
ICF. 

We conclude the following from our studies: 

1. T e (N + ) and T e (0 + ) diverge significantly from each other in low metallicity regimes. 
Therefore, the assumption of previous authors that these two temperatures are equiv- 
alent appears to be invalid at low metallicities. 

2. A realistic value for the N + /0 + ionization correction factor is 1.08±.09. 

3. Our final N/O values are systematically larger than previous results, while we find 
that O/H is generally below other measurements. This outcome is due to our new 
temperatures and ICF determinations. 

4. The majority (52) of objects in our survey form a relatively flat, horizontal system 
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in the N/O-O/H plane, consistent with the idea that both N and O production are 
primary at low metallicity. These galaxies comprise what we define as the N/O plateau. 

5. Plateau objects appear to form a Gaussian distribution in log(N/0) between -1.54 and 
-1.27 with the logarithm of the weighted arithmetic average, log(N/0), being -1.43 
(+.0084/-.0085) and a standard deviation of +.071/-.084. 

6. A x 2 analysis of the plateau objects indicates that only a small fraction of the observed 
scatter in N/ O is intrinsic, although this conclusion remains tentative until line strength 
uncertainties, which are the largest contributors to the abundance uncertainties, can 
be rigorously determined for a large sample of objects. 

7. The remaining 16 sample objects have log(N/0) values exceeding -1.27 and therefore 
reside above the N/O plateau as we define it here. 

Abundances of N and O are expected to improve in accuracy and precision as better 
temperature determinations become possible through the use of large telescopes ca- 
pable of obtaining high S/N spectral observations of the auroral lines of N + and + . 
Until that time, photoionization models of increased sophistication can go a long way 
toward reducing the uncertainty, as we have shown here. 

We are greatful to E. Baron, J. X. Prochaska, Schaerer D., T. Hoffmann, and G. J. 
Ferland for helpful discussions and comments. This research is supported by NSF grant 
AST 03-07118 and NASA grant NAG5-3505 to the University of Oklahoma. 
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Table I. ATOMIC DATA 



Ion Energy Levels Transition Rates Collision Strengths (at 10 4 K) 

N+ Nist Nist Lennon & Burke 1994 

0+ Nist Nist McLaughlin & Bell 1994 

0+ 2 Nist Nist Lennon & Burke 1994 

S+ Flower 1983 Ramsbottom et al. 1996 
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Table 2. PHOENIX MODELS 



Star 


Metal 


log(M) 




T <=// 


log(g) 


R/Rq 


log(L) 




abundances 51 


(M yr- 1 ) 


(km s- 1 ) 


(10 4 K) 


(cm s _1 ) 


(R q =6.95x10 10 cm) 


(erg s" 1 ) 


1 


solar/5 


-6.80 


2330 


4.265 


4.00 


10.48 


38.954 (-3.517/0.960) 


2 


solar/20 


-7.28 


1940 


4.260 


4.00 


10.48 


39.080 (-3.517/0.960) 


3 


solar /50 


-7.40 


1727 


4.260 


4.00 


10.48 


39.042 (-3.517/0.960) 


4 


solar / 5 


-6.41 


2550 


5.000 


4.00 


9.80 


39.560 (-2.040/0.960) 


5 


solar/20 


-6.89 


2130 


5.000 


4.00 


9.80 


39.500 (-3.040/0.960) 


6 


solar/50 


-7.21 


1894 


5.000 


4.00 


9.80 


39.300 (-2.040/0.960) 



a Exccpt those of N and C which in order to match typical values measured in low metallicity systems, were set 
such that log(N/0)=-1.46 and log(C/O)=-0.7. 
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Table 3. CLOUDY MODELS 



Star a 


12 + log(0/H) 


n(H) 
(cm- 3 ) 


r 

(pc) 


log(U) 


Geometry 


1 


8.0 


10 


2 


-1.37 


spherical 


1 


8.0 


10 


5 


-2.17 


thick shell 


1 


8.0 


100 


2 


-2.37 


thick shell 


1 


8.0 


100 


5 


-3.17 


plane-parallel 


2 


7.4 


10 


2 


-1.14 


spherical 


2 


7.4 


10 


5 


-1.94 


thick shell 


2 


7.4 


100 


2 


-2.14 


thick shell 


2 


7.4 


100 


5 


-2.94 


thick shell 


3 


7.0 


10 


2 


-1.16 


spherical 


3 


7.0 


10 


5 


-1.96 


thick shell 


3 


7.0 


100 


2 


-2.16 


thick shell 


3 


7.0 


100 


5 


-2.96 


plane-parallel 


4 


8.0 


10 


2 


-0.39 


spherical 


4 


8.0 


10 


5 


-1.19 


spherical 


4 


8.0 


100 


2 


-1.39 


thick shell 


4 


8.0 


100 


5 


-2.19 


thick shell 


5 


7.4 


10 


2 


-0.42 


spherical 


5 


7.4 


10 


5 


-1.22 


spherical 


5 


7.4 


100 


2 


-1.42 


thick shell 


5 


7.4 


100 


5 


-2.22 


thick shell 


6 


7.0 


10 


2 


-0.68 


spherical 


6 


7.0 


10 


5 


-1.48 


spherical 


6 


7.0 


100 


2 


-1.68 


thick shell 


6 


7.0 


100 


5 


-2.48 


thick shell 



a The nebular and stellar abundances are identical for all ele- 
ments except for Na, Mg, Al, Si, Ca and Fe, whose abundances 
are 0.1 X their stellar value. 
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Table 4. ERROR PROPAGATION 



Quantity Contributing Factors 



R[0 III] ^4363, ^4959 > ^5007 

T e (0+ 2 ) R[0 III] 

T e (0+) T e (0+ 2 ) 

T e (N+) T e (0 + 2 ) 

0+ 2 /H+ /5007, T e (0+ 2 ) 

O+/H+ / 3 727, T e (0+ 2 ) 

N+/H+ / 6584 , T e (0+ 2 ) 

O/H 0+, 0+ 2 

N/O ICF, N+/H+, O+/H+ 



Table 5. PHYSICAL CONDITIONS AND IONIC ABUNDANCES 



Object Names N e T e (0+ 2 ) T e (0+) T e (N+) O+/H+ 0+ 2 /H+ N+/H+ 

(cm- 3 ) (K) (K) (K) (xlO 5 ) (xlO 5 ) (xlO 5 ) 



T ZW IS NW MRK 116 0930+554 


<50 


19000 


+_ 


800 


15800 


+_ 


100 


14900 


+_ 


100 


o 


.178 


+_ 


0. 


008 


I 


338 




o 


.080 


o 


006 


+_ 


0. 


000 


I ZW 18 SE, MRK 116, 0930+554 


200 


16700 


+ 


1400 


15300 


± 


500 


14700 


+ 


300 


0. 


.365 


± 


0. 


.045 


1 


.539 


± 





.298 


o 


Oil 


+ 


0. 


001 


SRS 0335-052 


250 


20800 


+_ 


400 


15900 




100 


14600 


+_ 


100 


o 


.201 




0. 


003 


I 


.802 


± 


o 


.021 


o 


007 


+_ 


0. 


002 


SRS 0040+544 N 


150 


22100 


+ 


600 


15600 


± 


100 


14100 


+ 


300 


o 


.364 


± 


0. 


017 


1 


.971 


± 


o 


.028 


0. 


.014 


+ 


0. 


002 


SRS 1 1 594-545 


100 


19700 


+ 


600 


15900 


± 


100 


14900 


+ 


100 





.460 


± 


0. 


014 


2 


138 


± 


o 


.076 


o 


.017 


+ 


0. 


002 


SRS 1415+437 


50 


18600 


+ 


300 


15800 


± 


100 


14900 


+ 


100 





725 


± 


0. 


012 


1 


.996 


± 


o 


.047 


o 


.025 


+ 


0. 


001 


KTSSR 85 


750 


18100 


+_ 


1100 


15700 




200 


14900 


+_ 


100 


o 


819 


± 


0. 


063 


2 


.258 


± 


o 


.259 


o 


033 


+_ 


0. 


005 


TTCJC 4483 0832+699 


50 


16700 


+_ 


300 


15300 


± 


100 


14700 


+_ 


100 


o 


.803 




0. 


023 


2 


.321 




o 


.086 


o. 


023 


+_ 


0. 


001 


TOL 1214-277 


200 


19100 


+ 


700 


15800 


± 


100 


14900 


± 


100 


o 


.158 


± 


0. 


012 


3 


.265 


± 


o 


.178 


o 


.007 


+ 


0. 


001 


KISSR 1490 


200 


19300 


± 


1100 


15900 


± 


100 


14900 


± 


100 


1 


.046 


± 


0. 


049 


2 


.465 


± 





.210 


0. 


.048 


± 


0. 


.008 


SBS 1211+540 


150 


17200 


± 


200 


15400 


± 


100 


14800 


± 


100 





485 


± 


0. 


013 


3 


.605 


± 





.113 





016 


± 


0. 


001 


KISSR 1013, KISSR 211 


400 


18200 


± 


1100 


15700 


± 


200 


14900 


± 


100 


1 


.573 


± 


0. 


117 


2 


.709 


± 





.306 


0. 


.117 


± 


0. 


010 


SRS 1249+493 


200 


16900 


± 


700 


15300 


± 


200 


14700 


± 


200 





827 


± 


0. 


055 


3 


.517 


± 





.333 





.026 


± 


0. 


003 


116+583 R 


500 


16800 


± 


900 


15300 


± 


300 


14700 


± 


200 





.490 


± 


0. 


040 


3 


.922 


± 





.491 





.021 


± 


0. 


.006 


VII Zw 403, 1124+792 


<50 


15600 


± 


200 


14800 


± 


100 


14300 


± 


100 


1 


128 


± 


0. 


030 


3 


.397 


± 





.130 





.042 


± 


0. 


002 


TOL 1304-353 


200 


18600 


± 


700 


15800 


± 


100 


14900 


± 


100 





.258 


± 


0. 


009 


4 


.527 


± 





.373 





Oil 


± 


0. 


003 


SRS 1420+544, CG 413 


200 


18000 


± 


400 


15700 


± 


100 


14900 


± 


100 





.402 


± 


0. 


013 


4. 


419 


± 





.174 





.014 


± 


0. 


.002 


SRS 1205+557 


200 


16200 


± 


700 


15100 


± 


300 


14600 


± 


200 


1 


759 


± 


0. 


121 


3 


.312 


± 





.345 





.067 


± 


0. 


.005 


SRS 1128+573 


200 


17000 


± 


600 


15400 


± 


200 


14800 


± 


100 





.668 


± 


0. 


035 


4. 


587 


± 





.387 





.026 


± 


0. 


005 


C 1543+091 


<50 


16500 


± 


500 


15200 


± 


200 


14600 


± 


100 





.448 


± 


0. 


024 


4. 


.812 


± 





.403 


0. 


.026 


± 


0. 


006 


UM 461, SCHG 1148-020 


50 


17200 


± 


400 


15500 


± 


100 


14800 


± 


100 





.353 


± 


0. 


016 


4. 


.934 


± 





.290 





.012 


± 


0. 


.001 



Table 5 — Continued 



Object Names N e T e (0+ 2 ) T e (0+) T e (N+ ) O+/H+ 0+ 2 /H+ N+/H+ 

(cm- 3 ) (K) (K) (K) (xlO 5 ) (xlO 5 ) (xlO 5 ) 



I ZW 36, MRK 209, HARO 29, 1223+487 


50 


16300 


± 


100 


15100 


± 


100 


14600 


± 


100 


o 


.566 


+ 





.004 


4 


.865 


+ 





.048 





.023 


± 





001 


SBS 1331+493 N 


150 


16200 


± 


300 


15100 


± 


100 


14600 


± 


100 


o 


.723 


± 





.025 


4 


.862 


± 





.205 


0. 


026 


± 





002 


MRK 1434, SBS 1030+583 


200 


15600 


± 


100 


14800 


± 


100 


14300 


± 


100 


o 


.849 


± 





.016 


4 


.929 


± 





.121 





.026 


± 


0. 


001 


MRK 193, SBS 1152+579 


150 


16400 


± 


200 


15200 


± 


100 


14600 


± 


100 





.709 


± 





016 


5 


.204 


± 





145 





032 


± 





002 


MRK 71, NGC 2343, 0723+692 B 


200 


15000 


± 


300 


14600 


± 


200 


14100 


± 


200 


1 


.483 


± 





.071 


4 


.637 


± 





.285 





.048 


± 


0. 


003 


MRK 600, 0248+042 


50 


15900 


± 


200 


14900 


± 


100 


14500 


± 


100 


1 


.069 


± 





.025 


5 


.116 


± 





.168 





029 


± 





001 


MRK 71, NGC 2343, SBS 0749+568 B 


200 


15300 


± 


800 


14600 


± 


400 


14200 


± 


300 


1 


.507 


± 





162 


4 


.980 


± 





.711 





.065 


± 


0. 


007 


MRK 36, 1102+294 


550 


14900 


± 


400 


14400 


± 


200 


14000 


± 


200 


o 


.808 


± 





.058 


5 


.722 


± 


0. 


414 





031 


± 


0. 


.007 


POX 4: C1148-203 


<50 


15100 


± 


400 


14500 


± 


200 


14100 


± 


200 


o 


.755 


± 





.054 


5 


.800 


± 





.425 





.041 


± 





002 


IVlxxrv 1410, 3dU Uyl/+DZ/ 


onn 
zUU 


1 e^onn 
10ZUU 


i 

± 


Qnn 


140UU 


i 

± 


onn 
ZUU 


1 a onn 
14ZUU 


± 


i nn 
1UU 


1 

1 


7Q A 

■ 1 o4 


± 


n 

u. 


n70 
Viz 


A 

4 


qqc; 
.OOO 


i 

± 


n 

u. 


zoo 


n 

u. 


nc;n 


± 


n 

u 


nno 
UUZ 


MRK 71, NGC 2343, 0723+692 A 


100 


15900 


± 


100 


14900 


± 


100 


14500 


± 


100 





.505 


± 


0. 


.003 


6 


.132 


± 


0. 


.054 


0. 


.020 


± 


0. 


001 


SBS 1533+574 A 


<50 


14500 


± 


600 


14100 


± 


400 


13800 


± 


300 


2 


.431 


± 


0. 


226 


4 


.573 


± 


0. 


505 


0. 


112 


± 


0. 


007 


MRK 1486, SBS 1358+576 


<50 


14800 


± 


200 


14300 


± 


200 


14000 


± 


100 


1 


.575 


± 





.054 


5 


.530 


± 





.235 





.096 


± 


0. 


003 


SBS 1331+493 A 


200 


13600 


± 


900 


13500 


± 


700 


13300 


± 


600 


2 


.837 


± 


0. 


.506 


4 


.278 


± 





803 





.102 


± 


0. 


013 


KISSR 675, KISSB 187 


1000 


15200 


± 


1100 


14600 


± 


600 


14200 


± 


500 


1 


.356 


± 


0. 


218 


5 


.901 


± 


1 


211 


0. 


.048 


± 


0. 


013 


FAIRALL 30 


100 


14600 


± 


300 


14200 


± 


200 


13900 


± 


100 


1 


.079 


± 





.044 


6 


.418 


± 


0. 


.373 





.059 


± 





003 


C 0840+120 


200 


14200 


± 


400 


13900 


± 


300 


13700 


± 


200 


1 


.540 


± 


0. 


111 


6 


.080 


± 


0. 


482 


0. 


.064 


± 





009 


SBS 0926+606 


150 


14400 


± 


200 


14100 


± 


200 


13800 


± 


100 


1 


.851 


± 


0. 


.072 


5 


.812 


± 





.262 


0. 


.076 


± 


0. 


002 


I ZW 49, MRK 59, NGC 4861, ARP 209, 1256+351 


100 


14600 


± 


700 


14200 


± 


500 


13900 


± 


400 


2 


.280 


± 


0. 


278 


5 


.460 


± 


0. 


767 


0. 


108 


± 





010 


KISSR 1194 


50 


14600 


± 


400 


14200 


± 


200 


13900 


± 


200 


1 


.976 


± 





.175 


5 


.935 


± 





.487 


0. 


.078 


± 


0. 


.008 


KISSR 396, KISSB 145, WAS 81 


200 


14100 


± 


500 


13900 


± 


300 


13600 


± 


300 


2 


.705 


± 


0. 


232 


5 


.352 


± 


0. 


514 


0. 


.078 


± 





005 



Table 5 — Continued 



Object Names N e T e (0+ 2 ) T e (0+) T e (N+) O+/H+ 0+ 2 /H+ N+/H+ 

(cm- 3 ) (K) (K) (K) (xlO 5 ) (xlO 5 ) (xlO 5 ) 



UM 420 0218+003 


200 


13900 


+ 


800 


13800 


+ 


600 


13500 


+ 


500 


2 


.754 


± 





.425 


5 


.384 


+ 





.915 


0. 


.267 


+ 


0. 


026 


1437+370 


200 


14200 


± 


300 


13900 


± 


200 


13600 


± 


200 


1 


.511 


+ 


0. 


.079 


6 


699 


± 


o 


.382 


0. 


.056 


± 


0. 


003 


KISSR 1778 


50 


13100 


+ 


900 


13100 


± 


700 


13000 


± 


700 


3 


.622 


± 


o 


816 


4, 


605 


± 


o 


.968 


o 


.149 


± 


0. 


027 


UM 462, SCHG 1150-021 


200 


14300 


± 


200 


14000 


± 


100 


13700 


± 


100 


1 


.527 


± 


0. 


.050 


6 


841 


± 





.275 


0. 


055 


± 





007 


SBS 1222+614 


200 


14300 


± 


100 


14000 


± 


100 


13700 


± 


100 


1 


.231 


± 





.026 


7, 


.369 


± 





.181 





.035 


± 


0. 


001 


TOL 1304-386 


100 


14200 


± 


300 


13900 


± 


200 


13600 


± 


200 


1 


.073 


± 


0. 


061 


7. 


.781 


± 


o 


.502 


0. 


082 


± 





006 


1533+469 


<50 


13900 


± 


500 


13800 


± 


400 


13500 


± 


300 


2 


.456 


± 


o 


.257 


6 


.514 


± 


o 


.740 





.149 


± 


0. 


012 


SBS 0907+543 


100 


14500 


± 


400 


14200 


± 


300 


13900 


± 


200 





.939 


± 


0. 


069 


8. 


113 


± 


o 


689 


0. 


030 


± 





005 


UM 469 


200 


12900 


+ 


1200 


13000 


± 


1000 


12800 


± 


900 


2 


.458 


± 


o 


710 


6 


.662 


± 


1 


.858 





.201 


+ 





039 


1054+365, CG 798 


200 


13900 


± 


200 


13700 


± 


100 


13500 


± 


100 


1 


.210 


± 


0. 


.044 


8. 


.010 


± 





.315 


0. 


.047 


± 


0. 


.002 


MRK 1450, SBS 1135+581 


200 


13200 


± 


200 


13200 


± 


100 


13000 


± 


100 


1 


.764 


± 





.063 


7. 


.459 


± 





.258 





.071 


± 


0. 


002 


MRK 1271, TOL 1053+064 


50 


14200 


± 


200 


13900 


± 


100 


13600 


± 


100 


1 


.677 


± 


0. 


055 


7. 


752 


± 





.285 


0. 


.086 


± 





003 


NGC 5253-6 


100 


12700 


± 


600 


12800 


± 


500 


12600 


± 


500 


4 


.204 


± 


0. 


.662 


5. 


.374 


± 





.793 


0. 


.284 


± 


0. 


029 


MRK 1304, UM 448, SCHG 1139+006 


100 


12200 


± 


600 


12400 


± 


500 


12300 


± 


500 


4 


.559 


± 





.766 


5. 


.090 


± 





.756 





.496 


± 


0. 


050 


MRK 22, SBS 0946+558 


50 


13500 


± 


200 


13500 


± 


200 


13300 


± 


100 


1 


.753 


± 


0. 


078 


7. 


.896 


± 





.354 


0. 


.064 


± 


0. 


003 


1441+294, CG 1258 


200 


13200 


± 


800 


13200 


± 


600 


13000 


± 


500 


1. 


.969 


± 


0. 


.340 


7. 


.709 


± 


1 


.326 





.087 


± 





012 


KISS 49, KISSB 94, CG 177 


100 


13100 


± 


800 


13100 


± 


700 


13000 


± 


600 


4 


.490 


± 


0. 


888 


5. 


.229 


± 


1 


.010 


0. 


239 


± 





037 


MRK 1409, SBS 0741+535 


550 


13700 


± 


700 


13600 


± 


500 


13400 


± 


400 


3 


.673 


± 





.486 


6 


.132 


± 





.869 





.121 


± 


0. 


011 


0948+532 


100 


13500 


± 


300 


13400 


± 


200 


13200 


± 


200 


1 


.631 


± 


0. 


096 


8. 


.250 


± 





.492 


0. 


.069 


± 





.004 


TOL 1345-420 


100 


13400 


± 


700 


13400 


± 


500 


13200 


± 


400 


2 


.034 


± 


0. 


.279 


7. 


.939 


± 


1 


.149 





.055 


± 


0. 


008 


UM 439 


600 


14200 


± 


300 


14000 


± 


200 


13700 


± 


200 





.890 


± 


0. 


.047 


9. 


.533 


± 





.579 


0. 


034 


± 





006 



Table 5 — Continued 



Object Names 


N e 


T e (0+ 2 ) 


T e (0+) 


Te(N+) 


O+/H+ 


0+ 2 /H+ 


N+/H+ 




(cm- 3 ) 


(K) 


(K) 


(K) 


(xlO 5 ) 


(xlO 5 ) 


(xlO 5 ) 


TOL 0645-376 


350 


12900 ± 1600 


13000 ± 1300 


12800 ± 1200 


3.272 ± 1.222 


7.342 ± 2.632 


0.254 ± 0.065 


KISSR 1845, CG 903, HS 1440+4302 


50 


13300 ± 300 


13300 ± 200 


13100 ± 200 


2.997 + 0.286 


8.061 + 0.622 


0.111 + 0.011 


MRK 5, 0635+756 


200 


12200 ± 500 


12400 ± 500 


12200 ± 400 


3.566 ± 0.515 


7.512 + 0.951 


0.168 + 0.016 


II ZW 40 


200 


13300 ± 200 


13300 ± 100 


13100 ± 100 


0.581 ± 0.026 


10.780 + 0.523 


0.057 ± 0.003 


MRK 1089, 0459-043 


100 


11100 ± 700 


11300 ± 700 


11200 ± 700 


5.533 ± 1.432 


6.302 + 1.249 


0.519 + 0.079 



Table 6. ELEMENTAL ABUNDANCES 



Object Names Ref. for 12 + log(§) logCyj) 12 + log(§) log(§) Ref. for 

/(A) (present work) (present work) (literature) (literature) cols. 5 & 6 



I ZW 18 NW, MRK 116, 0930+554 

I ZW 18 SE, MRK 116, 0930+554 

SBS 0335-052 

SBS 0940+544 N 

SBS 1159+545 

SBS 1415+437 

KISSR 85 

UGC 4483, 0832+699 
TOL 1214-277 
KISSR 1490 
SBS 1211+540 
KISSR 1013, KISSB 211 
SBS 1249+493 
116+583 B 

VII Zw 403, 1124+792 

TOL 1304-353 

SBS 1420+544, CG 413 

SBS 1205+557 

SBS 1128+573 

C 1543+091 

UM 461, SCHG 1148-020 



4 


7101 4. n Q23 

1 .101 _l_ U.WiitJ 


_i 


41 f\ + 049 


4 


7 2S(1 + 069 


_i 


ko5 + n 070 


10 


7 qn2 + 005 


_1 


40« + n 1 04 


5 


7 368 + 006 


_1 


382 + 071 


5 


7 415 + 013 


_1 


397 + 055 


6 


7 435 + n 008 


_1 


423 + 039 


12 


7 aqq + n 0^8 


_i 


+ 078 


5 


7 405 4- n 01 2 




509 + 041 


1 3 


7 534 + 023 


_1 


324 + 069 


12 


7.546 ± 0.027 


-1 


.300 ± 0.087 


5 


7.612 + 0.012 


-1 


.458 ± 0.044 


12 


7.632 ± 0.033 


-1 


.095 ± 0.061 


6 


7.638 + 0.034 


-1 


.468 ± 0.071 


8 


7.645 ± 0.048 


-1 


.326 ± 0.124 


8 


7.656 + 0.013 


-1 


.394 ± 0.042 


1 


7.680 + 0.034 


-1 


.316 ± 0.122 


6 


7.683 + 0.016 


-1 


.429 ± 0.082 


8 


7.705 + 0.031 


-1 


.382 ± 0.056 


8 


7.721 + 0.032 


-1 


.377 + 0.090 


1 


7.721 + 0.033 


-1 


.208 ± 0.117 


1 


7.723 + 0.024 


-1 


.417 + 0.050 



7. 


l™±g;g|g 


-1 


,560+g 


090 


7 


7. 


26o±8:81S 


-1 


.600±g 


.060 
060 


7 


7 


29o±8:8iS 


-1 


580±g 


.030 


11 


7. 


oic_l_0.016 
(1 016 


-1 


.523±g 


.068 
OS 1 


7 


7, 


■442±°;gi| 


-1 


,533+g 


.062 
072 


7 


7, 


504±g.oi4 


-1 


.568±g 


042 
47 


7 


7, 


500±g;g60 


-1 


380±g 


.110 
.110 


12 


7, 


ceq+0.015 


-1 


629±g 


053 


7 


7. 


K 7fi 4-0.025 


-1 


,449+g 


.080 

(i')N 


7 


7, 


560+g-gZ , 


-1 


410+P, 


130 
130 


12 


7, 


«7l +0.010 


-1 


.601+° 


.058 
.067 


7 


7, 


DDU± 0.050 


-1 


,180+g 


.090 
.090 


12 


7, 


703±8:B1I 


-1. 


.584+0 


.082 
101 


7 


7, 


6S0±8;8I8 


-1 


,450+g 


110 
110 


11 


7, 


fi qn+o.oio 

Dau=t 0.010 


-1 


530±g 


.030 
.030 


11 


7. 


791+0.033 
,zl=c 0.036 


-1 


.442+0 


.127 
181 


7 


7. 


72«+0.014 
' zo=l= 0.014 


-1 


.542+0 


.078 
.096 


7 


7. 


7*0+0.030 
' ou=l= 0.030 


-1 


,500+g 


.080 
.080 


11 


7, 


7^0+0-030 
' ou=c 0.030 


-1 


.510±g 


.070 
.070 


11 


7, 


77Q+0.025 
' ' o=E: 0.027 


-1 


,340+g 


.115 
.156 


7 


7, 


7R4+0.019 
' D4=t 0.020 


-1 


,579+g 


.059 
008 


7 



Table 6 — Continued 



Object Names 


Ref. for 


12 + log(§) 


log(o) 


12 + log(^) 


Mtj) 


Ref. ; 




/(A) 


(present work) 


(present work) 


(literature) 


(literature) 


cols. 5 


T 7,W 36 MRK 209 HARO 29 1223+487 


8 


7 7qk -|_ n 904 


_i qao -|_ n 039 

J_. Ortiz _1_ L/ . \) i J C7 


7 770+0.010 


_i 400+0 010 


11 


SRS 1 331 +493 N 


5 


7 747 + n nifi 

1 . 1 rfc 1 _1_ U-U1U 


_i 41 n + n n48 


7 8m +0.011 


-1 492+0- 050 


7 


MRK 1434 SRS 1030+583 


8 


7 7«o -|_ n 009 

1 . 1 Ui _1_ KJ . KJ \ ) O 


_i 470 -|_ n 049 


7 700+O.OIO 
'• ,yUx: 0.010 


-1 600± 02 

I.OUUICq Q2Q 


1 1 


MRK 193 SRS 1152+579 


5 


7 772 + 011 


— 1 310 + 043 


7 865±S-S 1 S 


-1 274±9'PJZ 
- L -^'^ :c 0.052 


7 


MRK 71, NGC 2343, 0723+692 R 




7 787 + 021 


— 1 456 + 049 


7 810± 0.020 
' •°- LUX: 0.020 


-1 570±S-S4S 
J..JI u=c 0.040 


11 


MRK 600, 0248+042 


10 


7.791 ± 0.012 


— 1.535 ± 0.040 


7 830±R R}R 

' - ooux: 0.010 


-1 670±K-P.§P. 

i,ulux .030 


11 


MRK 71 NGC 2343 SRS 0749+568 R 


8 


7 812 + 049 


— 1 332 ± 075 


7 g 50± 0.050 
' - OJU ^0.050 


-1 440±0}}0 


11 


MRK 36, 1102+294 


1 


7.815 ± 0.028 


— 1.382 ± 0.111 


7 847±P,-R22 
' ■ cv±l ^0.023 


-1 459±S495 


7 


POX 4: C1148-203 


1 


7.817 ± 0.028 


— 1.233 ± 0.052 


7 855±S'S2§ 

' - OJJ ^0.025 


-1 350±SSIi 


7 


MRK 1416, SRC 0917+527 


8 


7.821 ± 0.018 


— 1.507 ± 0.044 


7 860±P.'828 

1 .OUUXQ Q2Q 


-1 620±°-°iP 


11 


MRK 71, NGC 2343, 0723+692 A 


8 


7 822 + 004 


— 1 373 ± 040 


7 050+O.OIO 


-1 520± - 010 


11 


SRS 1 533+574 A 


8 


7 845 + 034 


— 1 303 ± 061 


7 880± ' 030 
' - oou ^0.030 


_i 430+0.080 


11 


IVlIXlv l^toO, ODO 1000tO(D 


o 


7 qco _i_ n m ^ 


1 1 70 _|_ fi 0A9 


7 aao+0.010 
'• 8 ° U± 0.010 


1 Q1O+0-030 


1 1 
1 1 


SRS 1331+493 A 


6 


7.852 ± 0.058 


-1.408 ±0.101 


7 aqq+0.047 


-1.495±g;?f| 


7 


KISSR 675, KISSR 187 


12 


7.861 ±0.074 


-1.418 ±0.140 


7 070+O.O9O 
'• 8 < U± 0.090 


-1.390±§;il8 


12 


FAIRALL 30 


1 


7.875 ±0.022 


-1.228 ±0.046 


7 g n7 _|_0.020 
'■ au,=c 0.021 


-i.339±8:8?g 


7 


C 0840+120 


1 


7.882 ±0.028 


-1.347 ±0.077 


7 01 riO.025 
'• yioi 0.027 


-1-448±8;° 9 1 


7 


SBS 0926+606 


8 


7.884 ±0.015 


-1.349 ±0.042 


7 010+O OIO 

' - ylu=l: o.oio 


_l 47O+0.030 
L4IUx .030 


11 


I ZW 49, MRK 59, NGC 4861, ARP 209, 1256+351 


13 


7.889 ±0.046 


-1.290 ±0.075 


7 040+0. 030 
'■ y4u=t 0.030 


-1 440+0-050 


7 


KISSR 1194 


12 


7.898 ± 0.028 


-1.367 ±0.068 


7 090+0. 040 
'• yzui 0.040 


_i 4KO+0 80 


12 


KISSR 396, KISSR 145, WAS 81 


12 


7.906 ± 0.030 


-1.508 ±0.058 


7 qoo+0 040 
' - yzu=t 0.040 


_1 coo+0.070 


12 



Table 6 — Continued 



Object Names Ref. for 12 + log(§) log(g) 12 + log(§) log(g-) Ref. for 

/(A) (present work) (present work) (literature) (literature) cols. 5 & 6 



UM 420, 0218+003 
1437+370 
KISSR 1778 

UM 462, SCHG 1150-021 

SBS 1222+614 

TOL 1304-386 

1533+469 

SBS 0907+543 

UM 469 

1054+365, CG 798 
MRK 1450, SBS 1135+581 
MRK 1271, TOL 1053+064 
NGC 5253-6 

MRK 1304, UM 448, SCHG 1139+006 

MRK 22, SBS 0946+558 

1441+294, CG 1258 

KISS 49, KISSB 94, CG 177 

MRK 1409, SBS 0741+535 

0948+532 

TOL 1345-420 

UM 439 



10 


7.911 ± 0.054 


— 0. 


.979 ± 0. 


087 


5 


7.914 ± 0.021 


— 1 


.400 ± 


050 


12 


7.915 ± 0.067 


— 1 


.351 ± 0. 


131 


1 


7.923 ± 0.015 


— 1 


.412 ± 


071 


8 


7.935 ± 0.009 


-1 


.508 ± 0. 


.039 


1 


7.947 ± 0.025 


-1 


.082 ± 0. 


.055 


6 


7.953 ± 0.038 


— 1 


.183 ± 0. 


.068 


8 


7.957 ± 0.033 


-1 


.461 ± 


.083 


1,3 


7.960 ± 0.095 


-1 


.052 ± 0. 


.156 


8 


7.965 ± 0.015 


-1 


.380 ± 


044 


5 


7.965 ± 0.013 


-1 


.363 ± 0. 


.041 


10 


7.975 + 0.013 


-1 


.257 + 


.042 


2 


7.981 + 0.047 


-1 


.136 + 0. 


.089 


10 


7.984 ± 0.048 


-0 


.929 ± 0. 


.092 


5 


7.984 + 0.016 


-1 


.400 ± 


044 


8 


7.986 ± 0.061 


-1 


.318 + 0. 


.103 


12 


7.988 ± 0.060 


-1 


.240 + 


.115 


8 


7.991 ± 0.044 


-1 


.448 ± 0. 


.078 


5 


7.995 + 0.022 


-1 


.339 ± 


.050 


1 


7.999 + 0.051 


-1 


.530 + 0. 


.092 


1 


8.018 + 0.024 


-1 


.385 ± 


.085 



7.930±g;8|g 


-1 


•080+°;!!° 


11 


7.965+0;°}! 


-1 


■ 4 88±aolo 


7 


7.930±B;8|8 


-1 


.440+°; !!° 


12 


7.960±o : g!| 


-1 


c nfi _i_0.080 
0.098 


7 


7.950+0;°!° 


-1 


.610±°;0J0 


11 


7.982±0;°2° 


-1 


■isi±8:8§4 


7 


8.042±0;°27 


-1 


•354±°;0?| 


7 


8.010±0;°3° 


-1 


■590±°;0|0 


11 


7.982±0;°|§ 


-1 


.113±° : }25 


7 


7.970+S-222 


-1 


.480±SS3S 

0.030 


11 


o nn7 _i_0.013 

o.uu(_c 013 


-1 


410+0.051 
•^^o.oss 


7 


7 qqn+0 010 


-1 




11 


a qqk-1-0.041 
°' UUD=C 0.045 


-1 


919-1-0.078 
zlz=c 0.096 


7 


7.990+0;°!° 


-1 


mn _|_0.100 
uiU:t 0.010 


11 


8.026±0;°!| 


-1 


4 qK_i_0.052 
• 4yo=t 0.059 


7 


7.990±0;°|° 


-1 


4 ln _|_0.130 
* iu=c 0.130 


11 


8.000+0;°!° 


-1 


OQQlO.110 


12 


8.010+0;°!° 


-1 


K 4 Q_|_0.100 

• 04U=t o.ioo 


11 


8.032+°;°!! 


-1 


.442±°;053 


7 


8.033±0;°3f 


-1 


•618±8;??l 


7 


8.045+°;°!° 


-1 


.444±°;0f0 


7 



Table 6 — Continued 



Object Names Ref. for 12 + log(§) log(g) 12 + log(§) log(-g) Ref. for 

I(X) (present work) (present work) (literature) (literature) cols. 5 & 6 



TOL 0645-376 


1 


8.026 ± 0.119 


-1.076 + 0.200 


8-046+0.?! 


_i ioq+0.148 


7 


KISSR 1845, CG 903, HS 1440+4302 


12 


8.044 + 0.027 


-1.397 + 0.071 


8.050±°;°i° 


_l 470+0.080 
1,,lltt 0.080 


12 


MRK 5, 0635+756 


10 


8.044 ± 0.042 


-1.292 + 0.083 


8.040±°;°t° 


_l o fi n+0.100 
l.ODUlt 1Q0 


11 


II ZW 40 


1 


8.055 ± 0.020 


-0.974 ± 0.048 


8.046±°;°|1 


-1.042±°;°lf 


7 


MRK 1089, 0459-043 


10 


8.073 ± 0.070 


-0.993 + 0.135 


8.040+°;°!° 


_i nKn+0.140 


11 



REFERENCES.-(l)Campbell ct al. 1986;(2)Walsh & Roy 1989;(3)Pagel et al. 1992;(4)Skillman & Kennicutt 1993;(5)Izotov 
ct al. 1994;(6)Thuan ct al. 1995;(7)Kobulnicky & Skillman 1996;(8)Izotov ct al. 1997; (9)Izotov & Thuan 1998a; (lO)Izotov & Thuan 
1998b; (ll)Izotov ct al. 1999;(12)Melbourne et al. 2004;(13)Kobulnicky & Skillman 1998. 
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Fig. 1. — Modeled spectra generated with PHOENIX (solid lines) and with WM- 
Basic (dashed lines) for two O-dwarfs having O/H = solar/20. The top panel features a 
43 kK star, while the bottom panel features a 50 kK star. 
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< Te(0++) > < Te(0+) > 

Fig. 2.— Left: T e (0+ 2 ) from R[0 III] (vertically) vs. T e (0+ 2 ) from eq. [6] (horizontally) for 
each model. Right: Parameterized T e (0 + ) (vertically) vs. T e (0 + ) from eq. [6] (horizontally). 
For the parameterized T e (0 + ), we used eq. [9] for the nlO models (open circles) and eq. [11] 
for the nlOO models (closed circles). 
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Fig. 3. — Electron temperature parameterizations given by our models (solid lines), and by 
pag92, izo94, and sta96 (dotted lines). The solid lines labeled eq. [8] (top panel) and eq. [9] 
(bottom panel) are best fits to nlO models (open circles). The solid lines labeled eq. [10] 
(top panel) and eq. [11] (bottom panel) are best fits to nlOO models (closed circles). T e ff, 
r D , and O/H are indicated for the models located at the two T e (0 +2 ) extrema. 



25000 




le+19 2e+19 3e+19 4e+19 

depth (cm) 

Fig. 4. — Example of how T e , N + /N, and + /0 vary with depth as n(H) increases. The two 
models shown have T e ff = 43 kK, r Q = 2 pc, O/H = solar/50, and they differ only in their 
value of n(H). The solid and dashed lines in the top panel give the T e -structure for the nlO 
and nlOO models (respectively), while they correspond to N + /N and + /0 in the middle 
(nlO model) and bottom (nlOO model) panels. The dotted lines lead to T e at Max{N + 
for each model. 
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Fig. 5. — Diagram showing T e (0 + ) vs. T e (N + ), where the temperatures were obtainted 
using line flux ratios from our 24 models. Models with n(H) = 10 cm -3 are shown with 
open circles, while models with n(H) = 100 cm~ 3 are shown with filled circles. The diagonal 
corresponds to points such that T e (N + ) = T e (0 + ). 




log(N+/0+), method 3 




log(N+/0+), method 3 

Fig. 6. — Comparison of N + /0 + ratios derived for our sample objects. The top panel shows 
method 1 vs. method 3, while the bottom panel shows method 2 vs. method 3. 




log(N+/0+) [present work, Ne=10cm A -3] log(N+/0+) [pag92, Ne=10cm A -3] 




log(N+/0+) [izo94, Ne=10 cm A -3] log(N+/0+) [sta96, Ne=10 cm A -3] 



Fig. 7. — Comparison of log(N + /0 + ) values derived for our sample objects using various 
parameterizations. The top-left panel is a plot of values obtained using Eqs. [10] & [11] 
(vertical axis) vs. values obtained using Eqs. [8] & [9] (horizontal axis). In the remaining 
panels values obtained using Eqs. [8] & [9] (vertical axis) are plotted against values obtained 
using pag92, izo94, or sta96, adopting N e = 10 cm~ 3 . 
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Fig. 8.— Difference log[ N+/0+ ] - log[ (N+/0+) ] as a function of log[ (N+/0+) ] for our 
CLOUDY models, where log[ N + /0 + ] is based on line flux ratios while log[ (N + /0 + ) } is 
based on eq. [12]. 
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Fig. 9. — Present work vs. published 12 + log(0/H) values for our sample of 68 low 
metallicity systems. The literature values were taken from the sources indicated in the 
footnote of Table 6. The diagonal represents points of one to one correspondance. 
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literature 

Fig. 10.— Same as Fig. 9 but for log(N/0). 
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literature 

Fig. 11. — Same as Fig. 9 but for the uncertainty in log(N/0). 
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Fig. 12. — Log(N/0) versus 12 + log(0/H) as determined in this study for our sample 
objects. The dashed lines show the upper and lower limits for the N/O plateau as we define 
them (see the discussion of Fig. 13). 
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Fig. 13. — Observed distribution of values of log(N/0) as determined in this paper. The 
vertical dashed lines show the upper and lower limits for the N/O plateau. 



